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SUMMARY 


Motion  base  driving  functions  required  for  simulating  helicopter  and  Vertical/Short 
Take-Off  Landing  (V/STOL)  aircraft  were  investigated,  developed  and  evaluated  on  the 
NAVAIRDEVCEN  three  degree-of-freedom  motion  base. 

The  motion  base  drive  investigation  was  carried  out  in  three  steps.  The  first  step 
entailed  a survey  of  both  the  civilian  and  military  branches  of  the  flight  simulation 
community  for  the  purpose  of  determining  the  state-of-the-art  moving  base  drive  tech- 
niqu(!S.  The  second  step  was  the  simulation  design  where  these  drive  techniques  were 
adapted  to  the  NAVAIRDEVCEN  simulation  facilities  and  Research  and  Development 
(R&D)  mission.  Finally,  the  motion  base  drive  functions  were  implemented  in  the  simu- 
lation laboratory  and  e:q)erienced  pilots  were  used  to  evaluate  and  optimize  the  drive 
techniques. 

Results  of  the  survey  indicated  that  all  drive  techniques  have  the  same  underlying 
assumption:  pilots  are  capable  of  sensing  two  kinds  of  motion  - rotational  acceleration 
and  specific  force  (the  apparent  force  at  the  pilot's  station).  Because  of  severe  motion 
base  acceleration,  velocity  and  position  constraints,  the  aircraft  motion  caimot  be  simu- 
lated exactly.  Washout  filters  must  be  employed  in  order  to  limit  the  travel  of  the 
motion  base  hardware.  Traditionally,  therefore,  motion  base  drive  functions  have  con- 
sisted of  washout  circuits  of  one  type  or  another. 

In  recent  years,  attempts  to  improve  motion  base  drive  functions  have  been  in 
three  areas:  1)  the  type  of  washout  filter  used,  2)  the  axis  system  in  which  washout 
filters  are  applied,  and  3)  the  use  of  a technique  known  as  "tilt"  to  simulate  low  fre- 
quency specific  forces. 

A hybrid  simulation  was  designed  to  investigate  these  three  areas.  In  particular, 
four  washout  filter  t3rpes  (first  order,  third  order,  high  order  with  tilt  and  nonlinear 
adaptive)  were  applied  to  the  aircraft  rotational  degrees  of  freedom  and  specific  force 
vector  in  body  axes,  inertial  axes  and  in  conjunction  with  a linearized  transformation 
from  body  axes  to  inertial  axes. 

An  experiment  was  devised  where  washout  filter  parameter  settings  were  chosen 
and  the  relative  merit  of  each  motion  base  drive  algorithm  was  evaluated. 

The  results  of  the  experiment  show  that  the  axis  system  of  application  for  washout 
has  no  effect  on  the  accuracy  with  which  pilots  fly  the  aircraft.  Indeed,  the  axis  sys- 
tem transformations  may  be  linearized  by  applying  the  small  angle  approximation 
without  jeopardizing  the  quality  of  the  motion. 
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High  order  washout  filters,  with  or  without  tilt,  produced  anomalous  cues  which 
led  to  pilot  induced  oscillations.  These  anomalous  cues  could  not  be  compensated  for 
in  the  motion  drive  algorithms  because  no  lateral  or  longitudinal  linear  drive  mecha- 
nism exists  for  a three  degree-of-freedom  motion  base. 

Pilots  flew  the  aircraft  with  greater  accuracy  when  first  order  washout  filters 
were  used  then  they  did  for  any  other  washout  t3^e.  Nonlinear  adaptive  washout  filters 
were  preferred  by  the  pilots  and,  indeed,  they  did  produce  the  most  accurate  recovery 
of  aircraft  motion.  The  poorer  performance  of  the  pilots  when  using  nonlinear  adaptive 
may  have  been  due  to  a turbulence  model  used  in  the  experiment.  Flight  in  turbulent 
air  produces  motion  which  actually  hinders  the  pilot  when  flying  the  aircraft.  The  non- 
linear adaptive  filter  was  more  effective  than  the  other  washout  filter  types  in  repro- 
ducing the  turbulent  motion  and  may  have  suffered  for  it. 
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DIS  CUSSI  ON 


BACKGROUND 

Piloted  flight  simulators  are  used  extensively  for  both  research  and  pilot  training 
purposes.  Escalating  costs  of  operation  and  aircraft  procurement  can  be  cited  as 
primary  reasons. 

It  has  been  clearly  established  that  an  important  aspect  of  piloted  flight  simulation 
is  associated  with  the  generation  of  realistic  motion  cues,  especially  in  simulations  of 
aircraft  with  stability  and  control  problems  such  as  the  V/STOL. 

Today's  state-of-the-art  motion  systems  are  classified  as  either  moving  base 
systems  or  sustained  force  systems.  A motion  platform  or  motion  base  is  an  example 
of  the  former,  while  a G-seat  or  manned  centrifuge  Is  classified  the  latter.  This  re- 
search project  was  concerned  with  investigating  the  drive  functions  for  a moving  base 
system. 

Motion  simulation  via  moving  base  simulators  is  made  possible  because  the  human 
nervous  system  is  highly  adaptable.  Specifically,  certain  types  of  stimuli  can  evoke 
the  same  sensations  and  responses  even  when  the  stimulus  magnitude  varies.  As  a 
result  moving  base  simulators,  which  contain  activators  and  drive  linkages  with  severe 
position,  velocity  and  acceleration  limitations,  are  not  called  on  to  produce  one-to-one 
motion  simulation  of  the  actual  aircraft.  The  best  motion  is  a compromise  which  at- 
tempts to  provide  the  simulator  test  pilot  with  those  flight  sensations  and  motion  cues 
that  he  needs  while  staying  within  the  constraints  of  the  motion  drive  system.  However, 
adaptability  of  the  nervous  system  is  limited  and  care  must  be  taken  to  "time"  the 
motion  stimuli  lest  disorientation  and  resulting  nausea  reactions  occur. 

The  functions  which  transform  aircraft  motion  to  motion  base  actuator  commands 
are  commonly  referred  to  as  "motion  cue  drive  function/logic."  These  functions  em- 
body the  knowledge  of  aircraft  dynamics,  control  theory,  and  human  sensory  systems. 

OBJECTIVES 

The  objectives  for  this  exploratory  development  program  were  to  investigate  mo- 
tion base  driving  functions  and  to  expand  and  enhance  e^ertise  at  NAVAIRDEVCEN  in 
the  area  of  simulating  motion  sensations  of  manned  aircraft,  ^ecifically,  this  project 
was  directed  toward  the  investigation,  development  and  evaluation  of  motion  base  driving 
functions  required  for  simulating  helicopter  and  V/STOL  aircraft. 
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APPROACH 

It  was  first  planned  to  survey  both  the  civilian  and  militaiy  branches  of  the  flight 
simulation  community  for  the  state-of-the-art  moving  base  techniques. 

The  second  step  was  the  simulation  design  where  these  drive  techniques  were 
adapted  to  the  NAVAIRDEVCEN simulation  facilities  and  R&D  mission.  Included  in  the 
second  step  of  the  project  was  the  refurbishment  and  dynamic  optimization  of  the 
NAVAIRDEVCEN  motion  base  hardware.  After  the  motion  base  drive  functions  were 
implemented  in  the  simulation  laboratory,  experienced  pilots  were  used  to  evaluate 
and  optimize  the  drive  techniques. 
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MOTION  BASE  DRIVE  INVESTIGATION 


GENERAL  DISCUSSION 

A bibliography  of  research  articles  dealing  with  the  physics  and  physiology  of  mo- 
tion simulation  via  motion  base  hardware  was  compiled  and  is  included  at  the  end  of 
this  report.  Highlighting  the  survey  were  visits  to  NASA  Ames  Research  Center, 
Moffett  Field,  California,  and  the  Naval  Training  Equipment  Center  (NTEC),  Orlando, 
Florida. 

It  is  apparent  from  the  survey  that  a variety  of  techniques  exist  for  simulating 
aircraft  motion  on  a motion  base  platform.  However,  all  of  the  various  drive  fimctions 
have  the  same  underlying  assumption.  Pilots  are  capable  of  sensing  two  kinds  of 
motion: 

1.  Rotational  Acceleration 

2.  Specific  Force 

The  first  is  self  explanatory.  Specific  force  is  the  apparent  force  at  the  pilot's 
station.  It  is  the  quantity  we  can  feel  by  virtue  of  surface  pressures  on  our  bodies. 
Therefore,  specific  force  is  that  which  can  be  sensed  through  a I'elative  deformation 
of  different  parts  of  the  body.  Specifically,  it  is  the  sum  of  the  external  forces  divided 
by  the  vehicle  mass  less  the  gravitational  components,  i.e.,  it  is  the  lift,  drag,  thrust, 
and  Coriolis  accelerations  at  the  pilot's  station. 

Simulation  of  rotational  acceleration  and  specific  force  on  a motion  base  requires 
the  use  of  washout  circuits.  Commonly,  these  washout  circuits  receive  as  an  input 
the  computed  aircraft  motion  and  provide  a motion  onset  cue  to  the  pilot  in  the  simu- 
lator. Low  frequency  accelerations  are  washed  out  to  prevent  the  motion  base  plat- 
form from  exceeding  its  position  constraints. 

In  the  early  days  of  motion  simulation,  the  drive  functions  accepted  aircraft  ac- 
celerations resolved  in  a body  axis  coordinate  system  and  passed  them  through  linear 
first  order  washout  filters.  The  washed-out  accelerations  were  then  transformed  to 
an  inertial  axis  system  where  they  were  doubly  integrated  and  converted  to  motion 
base  position  commands  (see  figure  1).  The  frequency  breakpoints  of  the  washout 
filters  were  adjusted  such  that  the  washout  rate  was  less  than  the  sensoi'y  threshold 
of  the  pilot.  As  a result,  the  pilot  ideally  felt  the  onset  of  the  motion  and  then  the 
platform  was  gi’adually  slowed  to  a halt  at  a rate  imperceptable  to  him.  Gains  were 
adjusted  to  insure  that  the  simulator  did  not  exceed  its  position  constraints. 
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FIGURE  1 - Basic  Washout  Circuit 
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Attempted  Improvements  to  the  original  motion  drive  functions  have  been  in  three 
areas;  1.)  the  choice  of  axis  system  in  which  the  washouts  are  performed,  2.)  the 
type  of  washout  filter  used  (first  order  linear,  high  order  linear,  nonlinear  adaptive, 
etc.),  and  3.)  the  use  of  "tilt"  to  simulate  low  frequency  specific  forces. 

Since  motion  is  sensed  in  a body  centered  frame  of  reference  it  is  aesthetically 
desirable  to  washout  the  accelerations  in  that  axis  system.  However,  anomalies  as- 
sociated with  washout  in  one  component  of  the  body  axis  accelerations  would  then  be 
coupled  into  three  components  when  they  are  transformed  into  the  inertial  frame  of 
reference.  This  coupling  can  produce  a degree  of  simulator  drift  during  operation. 
Washing  out  the  accelerations  after  they  have  been  transformed  to  inertial  axes  would 
eliminate  the  drift  problem.  However,  washing  out  accelerations  in  a frame  of  ref- 
erence other  than  the  axis  system  in  which  the  accelerations  are  sensed  may  have  an 
effect  on  the  quality  of  the  motion. 

Both  techniques  have  been  suggested.  It  has  also  been  suggested  that  the  small 
angle  approximation  be  employed  to  linearize  the  transformation,  thereby  obviating 
the  question  of  whether  to  washout  in  inertial  or  body  axes. 

The  second  area  of  modification  to  the  original  washout  circuits  was  in  the  design 
of  the  washout  filters  themselves.  First  order  linear  filters  have  the  advantage  of 
simplicity.  However,  they  also  have  some  disadvantages.  The  frequency  cutoff  is  not 
sharp  for  a first  order  filter.  Therefore  accelerations  at  frequencies  near  the  filter 
breakpoint  will  suffer  large  distortions. ‘ Another  disadvantage  is  associated  with  the 
steady  state  response  of  a first  order  filter  to  step  inputs.  If  an  acceleration  step  is 
input  to  the  filter  and  then  doubly  integrated,  the  steady  state  output  will  be  some 
nonzero  position.  For,  although  the  acceleration  is  washed  out,  the  position  is  not. 

A series  of  acceleration  steps,  therefore,  will  eventually  cause  the  simulator  to  e.\ceed 
its  position  constraints. 

Second  order  linear  filters  have  sharper  cutoffs  than  first  order  filters,  but  they 
also  will  command  a steady  state  nonzero  position  for  an  acceleration  step. 

Third  order  washout  filters  have  sharp  cutoffs  and  they  also  insure  that  the  posi- 
tion, as  well  as  the  acceleration,  will  return  to  zero.  However,  high  order  linear 
filters  can  have  undesirable  transient  responses. 

Other  washout  filters  in  use  today  are  of  the  nonlinear  adaptive  variety.  These 
nonlinear  filters  attempt  to  minimize  a cost  function  designed  to  produce  one-to-one 
motion  whenever  possible  and  to  scale  down  the  filter  gains  whenever  the  input  is  too 
large  for  the  constraints  of  the  motion  base.  Such  filters  can  produce  one-to-one 
motion  simulation  a large  percentage  of  the  time.  However,  it  may  be  argued  that 
linear  time  invariant  washouts  are  more  desirable  because  they  allow  the  pilot  to  judge 
different  cases  by  comparing  relative  changes  in  the  cues. 
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Rotation  angles  (tilt)  of  the  cockpit  alter  the  specific  force  on  the  pilot  since  the 
gravity  vector  is  rotating  in  a body  centered  frame  of  reference.  This  effect  is  often 
used  to  simulate  low  frequency  specific  forces  in  the  lateral  and  longitudinal  directions. 
Specific  forces  arising  from  tilt  can  be  held  indefinitely,  unlike  the  force  sensations 
generated  by  translational  drive  mechanisms.  The  common  technique  for  utilizing  tilt 
is  shown  in  figure  2. 

Small  rotation  rate  signals  are  generated  that  drive  the  simulator  toward  an  atti- 
tude which  gives  appropriate  force  cues.  The  rotation  rates  must  remain  below  the 
sensory  threshold  of  the  pilot  to  prevent  anomalous  rotation  cues.  It  is  desired  to 
achieve  a specific  force  sensation  by  placing  the  gravity  vector  gc  in  a direction  along 
the  desired  specific  force  V(!ctor  f^  . This  is  done  by  rotation  the  simulator  about  an 
axis  perpendicular  to  both  f^  and  g^.  until  they  are  aligned.  Figure  3 shows  a tyiiical 
circuit  for  implementing  tilt.  The  specific  force  vector  is  resolved  into  inertial  axes 
and  the  cross  product  fcX(-gc)  generates  the  rotation  vector  which  will  rotate  the 
simulator  to  null  the  cross  product  and  align  the  desired  specific  force  vector  with  the 
gravity  vector.  Note  the  provision  for  a shaping  filter  (a  low  pass  filter  of  some  sort) 
which  passes  only  slowly  varying  force  terms  to  prevent  the  generation  of  large  anom- 
alous rotation  rates.  The  tilt  rates  are  then  summed  with  the  rates  from  the  rotational 
channels.  For  a six  degree-of-freedom  simulator,  it  is  possible  to  coordinate  the 
translational  and  rotational  channels  and  reproduce  longitudinal  and  lateral  specific 
forces  at  all  frequencies.  The  translational  drive  reproduces  the  high  frequency  forces 
and  the  rotational  drive  takes  over  for  the  low  frequency  forces  through  the  use  of  tilt. 

SIMULATION  DESIGN 

A hybrid  simulation  design  was  selected.  The  digital  portion  of  the  progTam,  im- 
plemented on  a CDC  6600,  was  responsible  for  simulation  control  (real-time  operation 
and  input/output  functions),  and  the  motion  base  drive  logic.  The  aircraft  simulation 
resided  on  an  EAI  8800  hybrid  analog  computer  which  also  acted  as  an  interface  be- 
tween the  digital  motion  base  drive  logic  and  the  Singer- Link  three  degree-of-freedom 
motion  base.  A digital  turbulence  model  (detailed  in  Appendix  A)  was  added  to  increase 
the  realism  of  the  simulation.  Figure  4 depicts  a block  diagram  of  the  simulation  ele- 
ments and  interfaces. 

The  motion  base  with  its  helicopter  controls  and  displays  in  an  enclosed  cockpit 
and  the  aircraft  model,  a UH-1  helicopter,  were  borrowed  from  the  Light  Airborne 
Multi-P>urpose  System  (LAMPS)  simulation  In  residence  at  the  NAVAIRDEVCEN  simu- 
lation facility  (sec  figure  5). 

The  motion  base  hydraulics  were  refurbished  with  new  hydraulic  fluid  and  filters. 
The  accumulator  was  recharged  with  nitrogen  and  the  hydraulic  system  pressure  raised 
to  the  desired  oj)erating  level.  New  servo  valves  were  purchased  when  an  old  one  was 
found  defective,  and  the  electronic  gains  were  adjusted  to  obtain  optimum  response 
from  the  system.  Appendix  B describes  the  refurbishment  of  the  motion  system. 
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FIGURE  2 - Tilt  Geometry 
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FIGURE  4 - Simulation  Block  Diagram 
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FIGURE  5 - Motion  Base  With  Enclosed  Cockpit 
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The  motion  base  drive  logic  was  designed  such  that  washouts  could  be  performed 
in  inertial  axes,  body  axes  or  in  conjunction  with  a linear  transformation.  Washout 
types  were  also  selectable:  first  order  linear,  third  order  linear,  high  order  linear 
with  tilt  or  adaptive  nonlinear  filters. 

Because  the  motion  base  has  only  three  degrees  of  freedom  (pitch,  roll  aj\d  heave) , 
the  utility  of  tilt  could  not  be  fully  tested.  It  has  been  shown  in  the  past,  that  tilt  algo- 
rithms are  most  effective  when  linear  and  rotational  drives  are  "coordinated".  Rota- 
tional washouts  produce  anomalous  forces  which  can  be  eliminated  by  a coordinated 
linear  drive.  Also,  specific  forces  can  be  simulated  at  all  frequencies  by  allowing  the 
linear  drives  to  pass  high  frequency  forces  and  angular  drives  (through  tilt)  to  simulate 
low  frequency  forces. 

Beeause  the  NAVAIRDEVCEN  motion  base  does  not  have  a longitudinal  or  lateral 
linear  drive,  coordinated  drive  logic  could  not  be  examined.  Therefore,  the  only  ques- 
tion concerning  tilt  that  was  addressed  by  this  research  has  to  do  with  the  effectiveness 
of  using  tilt  in  the  absence  of  linear  drives. 

The  tilt  algorithm  used  in  this  study  was  derived  from  the  drive  logic  developed  by 
J.  B.  Sinacori  for  Northrop  Corporation  in  1973.  ^ 

The  nonlinear  adaptive  washout  filter  was  patterned  after  the  algorithm  developed 
by  Parrish,  Dieudonne,  Bowles  and  Martin  for  NASA  Langley  Research  Center  ui  1973.2 
It  applies  the  method  of  continuous  steepest  descent  to  define  adaptive  washout  param- 
eters such  that  a given  cost  function  is  optimized. 

EXPERIMENTAL  DESIGN 

A total  of  thirteen  algorithms  were  tested.  Each  of  the  four  washout  filter  types 
(first  order  linear,  third  order  linear,  high  order  with  tilt  and  nonlinear  adaptive) 
were  applied  in  body  axes,  inertial  axes  and  in  conjunction  with  linear  transformations 
making  a total  of  twelve  motion  algorithms  and  no  motion  was  the  thirteenth  algorithm. 

A description  of  each  of  the  twelve  motion  algorithms  is  included  in  Appendix  C. 

The  experiment  was  performed  in  three  phases.  In  the  first  phase  a NAVAIR- 
DEVCEN helicopter  pilot,  LCDR  A.  Bingley,  flew  the  simulator  without  motion  to  check 
out  the  helicopter  model  and  suggest  changes  to  make  the  model  more  realistic.  .After 

^Sinacori,  J.  B. ; "A  Practical  Approach  to  Motion  Simulation";  AIAA  Paper  No.  73-931; 

Sep  1973. 

2 

Parrish,  R.  V. , Dieudonne,  J.  E.,  Bowles,  R.  L. , and  Martin,  D.  J. ; "Coordinated 

Adaptive  Washout  for  Motion  Simulators";  AIAA  Paper  No.  73-930;  Sep  1973. 
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l.CDH  lUnuley  was  satisfied  with  the  helicopter  model,  he  flew  each  of  the  twelve 
motion  algorithms  for  the  purpose  of  setting  washout  parameters.  These  settings  were 
then  used  for  the  rest  of  the  experiment. 

The  second  phase  of  the  experiment  tested  each  of  the  thirteen  algorithms  - twelve 
different  motion  algorithms  and  one  with  no  motion.  Four  NAVAIRDEVCEN  helicopter 
pilots  were  used  in  the  second  phase  of  the  experiment:  LCDR  A.  Bingley,  LCDR  G. 
Montgomery,  LCDR  C.  Park  and  LT  F.  Porter. 

Th(-  pilrjts  were  given  a flying  task  designed  to  exercnsc-  ;dl  the  degrees  of  freedom 
of  the  aircraft.  Referring  to  figure  6,  the  task  was  a race  track  pattem  with  a climb 
and  descent.  The  run  was  started  with  the  aircraft  heading  north  at  ;ui  altitude  of  500 
feet  and  an  airspeed  of  120  knots.  The  pilot  was  tlien  asked  to  execute  a l!S0  degree 
turn  to  the  left  at  constant  altitude.  After  leveling  out  on  a southerly  heading,  he  was 
asked  to  climb  at  500  feet  per  minute  to  an  altitude  of  750  feet.  lie  then  leveled  out  and 
initiated  another  constant  altitude  180  degrees  turn  to  the  left.  After  leveling  out  on  a 
northerly  heading,  he  started  a descent  at  500  feet  per  minute  to  an  altitude  of  500  feet 
where  he  leveled  out,  completing  the  run.  Airspeed  was  to  remain  constant  throughout 
the  run  which  lasted  approximately  3-1/2  minutes. 

Steering  commands  for  flying  the  pattem  were  provided  by  an  instrument  mounted 
above  the  instrument  panel.  Ideally,  motion  algorithms  should  be  tested  using  a Visual 
Flight  Rules  (VFR)  rather  than  an  Instrument  Flight  Rules  (IFR)  task.  However,  no 
out-the-window  visuals  were  available  for  this  simulation.  The  steering  command  in- 
strument, therefore,  was  positioned  above  the  instrument  panel  requiring  the  pilot  to 
set  up  a scan  pattem  similar  to  one  he  would  use  when  flying  VFR. 

Referring  to  figure  7,  the  steering  command  instrument  consisted  of  a horizontal 
bar,  a vertical  bar  and  a ball.  The  horizontal  bar  was  driven  with  airspeed  error, 
the  vertical  bar  with  heading  error  and  the  ball  with  altitude  error. 

The  pilots  practiced  with  no  motion  until  they  were  familiar  with  the  instmment, 
the  flying  pattem,  and  the  aircraft.  Then  they  practiced  with  each  motion  algorithm 
until  they  were  familiar  with  it,  and  data  was  collected. 

The  motion  algorithms  were  numbered  as  follows: 

Algorithm  Number  Algorithm  Type 

1 First  Order  Washout  (Linear  Transformation) 

2 Third  Order  Washout  (Linear  Transfomiation) 

3 Nonlinear  Adaptive  Washout  (Linear  Transformation) 

4 High  Order  Washout  with  Tilt  (Linear  Transformation) 

5 First  Order  Washout  (Body  Axis) 
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FIGURE  6 - Racetrack  Pattern 
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FIGURE  7 - In  strum  ont  ProvidinK  Flight  Task  Cues 


NADC-77306-20 


Algorithm  Number  Algorithm  Type 

6 Thii’d  Order  Washout  (Body  Axis) 

7 Nonlinear  Adaptive  Washout  (Body  Axis) 

8 High  Order  Washout  with  Tilt  (Body  Axis) 

9 First  Order  Washout  (Inertial  Axis) 

10  Third  Order  Washout  (Inertial  Axis) 

11  Nonlinear  Adaptive  Washout  (Inertial  Axis) 

12  High  Order  Washout  with  Tilt  (Inertial  Axis) 

13  No  Motion 

The  order  of  appearance  of  these  algorithjns  during  the  e^qjeriment  was  r.uidomly 
selected  for  each  pilot  to  minimize  the  effect  of  a "learning  curve": 


Algorithm  No. 

Pilot  #1 

2 

Pilot  #2 

8 

Pilot  #3 

13 

Pilot  #4 

12 

Algorithm  No. 

13 

2 

3 

1 

Algorithm  No. 

5 

11 

4 

2 

Algorithm  No. 

7 

10 

8 

13 

Algorithm  No. 

8 

7 

5 

7 

Algoritlim  No. 

10 

6 

12 

11 

Algorithm  No. 

4 

5 

7 

9 

Algorithm  No. 

11 

13 

9 

5 

Algorithm  No. 

3 

3 

6 

G 

Algorithm  No. 

12 

12 

2 

8 

Algorithm  No. 

1 

9 

10 

4 

Algorithm  No. 

6 

1 

11 

3 

Algorithm  No. 

9 

4 

1 

10 

The  data  collected  were  the  means  and  variances,  computer  on-line,  for  each  of 
nine  error  parameters.  The  following  formulas  were  used: 

T 

error  (t)dt 


T 

(error  (t))^dt-(Mean)^. 


Mean  of  error 


■ i / 


Variance  of  error  = 


^ I 
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The  nine  error  parameters  were:  Three  inertial  position  errors,  three  attitude 
angular  errors  and  three  stick  positions  (x  and  y cylic  control  and  collective). 

The  measures  of  effectiveness  were  taken  to  be  the  variances  of  the  error  param- 
eters rather  than  the  means,  since  the  variances  of  the  errors  carry  more  information 
on  the  accuracy  with  which  the  pilot  is  flying. 

The  attitude  and  position  error  variances  are  a measure  of  total  system  performance. 
The  variance  of  the  control  positions  gives  insight  into  the  reason  for  differences  in  sys- 
tem performance. 

Preliminary  results  of  the  second  phase  of  the  ejqjeriment  indicated  that  the  error 
variances  were  nearly  equal  when  washouts  were  performed  in  inertial  axes,  body  axes, 
or  in  conjunction  with  linear  transformations.  First  order  linear  filters,  high  order 
linear  filters  with  tilt  and  adaptive  nonlinear  filters  produced  error  variances  that  were 
approximately  equal  while  substantially  higher  variances  were  observed  for  third  order 
linear  filters. 

On  the  basis  of  these  results,  a third  phase  of  the  experiment  was  designed.  It 
was  decided  to  perform  the  washouts  in  inertial  axes  to  eliminate  the  problem  of  drift. 
First  order  filters,  high  order  filters  with  tilt  and  adaptive  nonlinear  filters  were  chosen 
to  be  tested  in  this  phase  of  the  experiment. 

Three  NAVAIRDEVCEN  pilots  were  used  in  the  third  phase:  LCDR  B.  Bowen, 

LCDR  S.  Hoxie  and  LT.  F.  Porter,  ft  was  felt  that  a larger  sample  size  was  needed 
to  distinguish  the  differences  between  the  three  algorithms.  Therefore,  the  pilots  flew 
the  three  algorithms  five  times  each  in  the  following  order: 

Pilot  *1 

Nonlinear  Adaptive  Washout 

First  Order  Washout 

High  Order  Washout  with  Tilt 

Pilot  #2 

First  Order  Washout 

High  Order  Washout  with  Tilt 

Nonlinear  Adaptive  Washout 

Pilot  »3 

High  Order  Washout  with  Tilt 

Nonlinear  Adaptive  Washout 

First  Order  Washout 
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In  addition  to  the  nine  error  parameters  of  the  second  phase  of  the  experiment, 
data  was  also  collected  on  the  means  and  variances  of  three  angular  rate  errors. 

A fast  fourier  transform  was  performed  on  the  control  position  time  histories 
showing  the  frequency  content  of  the  stick  movements.  It  was  felt  that  differences  in 
the  frequency  content  could  aid  in  understanding  any  observed  variations  in  system  per- 
formance for  the  three  algorithms. 


23 


NADC-77306-20 


DISCUSSION  OF  RESULTS 


PHASE  I 

The  washout  parameters  selected  in  this  phase  of  the  experiment  and  used  in  the 
second  and  third  phases  are  listed  in  the  algorithm  descriptions,  Appendix  C. 

PHASE  II 

The  second  phase  of  the  experiment  investigated  washout  types  and  the  axis  system 
in  which  washouts  are  applied. 

The  three  control  variances  (x  and  y cyclic  and  collective)  were  added  to  form  a 
control  variance;  the  three  Euler  angle  error  variances  were  added  to  form  an  attitude 
error  variance;  and  the  x,  y and  z error  variances  were  added  to  form  position  error 
variance.  The  control,  attitude  and  position  error  variances  were  then  averaged  over 
the  four  pilots  for  washout  types  and  axis  system  of  application.  Figures  8 through  10 
show  the  average  control,  attitude  and  position  error  variances  for  the  four  different 
washout  types.  Figures  11  through  13  are  the  average  control,  attitude  and  position 
error  variances  for  the  three  choices  for  axis  system  of  application. 

Control,  attitude  and  position  error  variances  are  larger  for  third  order  washout 
filters  than  for  the  other  three  washout  filters  - first  order,  high  order  with  tilt  and 
nonlinear  adaptive.  This  difference  is  seen  to  be  more  pronounced  for  attitude  :md 
position  error  variances  than  it  is  for  control  variances,  although  the  trend  is  the 
same  for  all  three  measures  of  effectiveness.  The  levels  of  sig-nificance  for  this  re- 
sult were  .30  for  control  variances,  .08  for  attitude  error  variances  and  . 10  for  posi- 
tion error  variances. 

The  control,  attitude  and  position  error  variances  for  axis  system  of  application 
(figures  11  through  13)  show  no  consistent  trend.  The  level  of  significance  for  the  con- 
trol variance  difference  is  .39.  For  attitude  error  variance  difference  the  level  of 
significance  is  .76;  and  for  position  error  variance  difference,  the  level  of  significance 
is  .60. 

Figures  14  through  16  show  the  control,  attitude  and  position  error  variances  for 
runs  with  motion  versus  runs  with  no  motion.  Smaller  variances  are  evident  for  runs 
with  motion  for  control  and  attitude  errors,  while  that  trend  is  reversed  for  position 
error  variances.  Th(!  one-tailed  test  levels  of  significance  for  these  results  were: 

Control  Variance  - .24 

Attitude  Error  Variance  - . 39 

Position  Error  Variance  - .45 
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FIGURE  9 - Attitude  Error  Variance  For  Washout  Types 
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FIGURE  10  - Position  Error  Variance  For  Washout  Types 
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FIGURE  12  - Attitude  Error  Variance  For  Axis  System  Of  Application 
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FIGURE  13  - Position  Error  Variance  For  Axis  System  Of  Application 
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FIGURE  14  - Control  Variance  For  Motion  Versus  No  Motion 
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FIGURE  16  - Position  Error  Variance  For  Motion  Versus  No  Motion 
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Thi'  pilots  could  sense  no  difference  in  the  motion  cues  when  the  axis  system  of 
application  for  washouts  was  varied.  However,  they  felt  very  definite  differences  in 
the  cues  when  washout  filter  types  were  varied. 

Pilots  unanimously  preferred  the  nonlinear  adaptive  washout  filter  because  the 
turbulence  was  most  apparent  for  this  type  of  filter.  Turbulence,  they  felt,  was  the 
most  important  part  of  the  motion  for  increasing  the  realism  of  the  simulation.  The 
high  order  washout  with  tilt  was  regarded  as  being  almost  as  realistic  as  the  nonlinear 
adaptive  washout,  but  less  turbulence  was  evident.  The  first  order  washout  was  de- 
scribed as  being  "smooth",  indicating  a lack  of  turbulence,  but  the  pitch  and  roll  cues 
were  excellent.  Most  pilots  detected  a tendency  for  Pilot  Induces  Oscillations  (lUO) 
when  the  third  order  washout  filters  were  used. 

The  pilots  felt  that  motion,  in  general,  added  a great  deal  of  n alism  to  the  simu- 
lation. However,  they  felt  that  they  did  not  consciously  use  the  motion  cues  to  help 

them  with  the  flying  task. 

i 

Figures  C-13  thru  C-24  in  Appendix  C show  the  response  of  each  of  the  twelve 
motion  algorithms  to  a left  turn  and  a climb. 

The  axis  system  of  application  makes  no  apparent  difference  in  the  time  response 
of  the  four  washout  filter  t3rpes.  A significant  difference  is  seen,  however,  for  the 
time  response  of  each  of  the  washout  filter  types  viewed  separately. 

First  order  washout  filters  show  very  good  motion  recovery  for  the  angular  degrees 
of  freedom.  Simulator  longitudinal  specific  force  follows  the  aircraft  motion  marginally 
well,  but,  the  simulator  lateral  specific  force  seems  to  bear  no  relation  to  the  air- 
craft motion.  Vertical  specific  force  recovery  is  fair  but  necessarily  of  very  low 
amplitude  since  the  simulator  total  travel  in  the  vertical  direction  is  only  ±1  foot. 

For  third  order  washout  filters,  the  simulator  angular  motion  follows  the  general 
motion  of  the  aircraft,  but  the  angular  washout  characteristic  of  high  order  filters  re- 
sults in  the  simulator  body  rate  being  of  opposite  sign  to  the  aircraft  body  rate  on  oc- 
casion, thus  producing  erroneous  cues.  Simulator  longitudinal  specific  force  once  again 
follows  the  aircraft  motion  marginally  well,  while  the  simulator  lateral  specific  force 
does  not.  The  vertical  specific  force  is  produced  in  the  same  manner  as  for  first  order 
washouts. 

High  order  washout  filters  with  tilt  also  follow  the  general  aircraft  motion  for 
angular  degrees  of  freedom,  but  the  sign  reversal  evident  for  third  order  washouts 
still  occurs  when  tilt  is  added. 

The  longitudinal  specific  force  recovery  appears  to  be  not  much  better  than  for 
the  other  filters  without  tilt.  Lateral  specific  force  recovery  is  only  barely  recog- 
nizable. Apparently,  most  of  the  aircraft  specific  force  is  in  the  high  frequency  range 
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and  is  thus  being  filtered  out  by  the  tilt  shaping  filter.  Once  again  the  vertical  specific 
force  is  handled  in  the  same  manner  as  for  the  other  linear  washout  circuits. 

Angular  rate  recovery  is  excellent  for  nonlinear  adaptive  washout  circuits.  High 
frequency  rates  are  reproduced  exactly  with  washout  occurring  only  for  sustained  air- 
craft rates.  No  evidence  of  sign  reversal  is  apparent.  Simulator  longitudinal  specific 
force  follows  aircraft  motion  as  well  or  better  than  the  linear  washout  circuits.  Lateral 
specific  force  is  reproduced  much  better  with  the  nonlinear  adaptive  washouts  than  for 
any  of  the  linear  washouts  including  the  high  order  washout  with  tilt.  Vertical  specific 
force  recovery  is  very  poor.  For  this  degree-of- freedom,  the  linear  washouts  per- 
formed much  better. 

PHASE  III 

The  third  phase  of  the  experiment  further  investigated  the  three  washout  filter  typos 
that  appeared  nearly  equal  from  the  second  phase  - nonlinear  adaptive  washout,  high 
order  washout  with  tilt  and  first  order  washout. 

Control,  attitude  and  position  error  variances  were  formed  as  in  the  second  phase 
of  the  e;q)erlment.  In  addition,  the  body  axis  x,  y and  z angular  rate  error  variances 
were  added  to  form  a body  rate  error  variance. 

Figures  17  through  20  show  the  control,  body  rate,  attitude  and  position  error 
variances  averaged  over  all  the  pilots'  runs  for  each  of  the  three  washout  filter  types. 

A consistent  trend  is  seen  for  control,  body  rate  and  attitude  error  variances.  First 
order  washouts  show  the  least  error  variances,  high  order  washout  with  tilt  show  higher 
error  variances.  For  position  error  variance,  first  order  washout  is  still  best,  however 
nonlinear  adaptive  washouts  show  better  position  error  variance  than  the  high  order 
washout  with  tilt.  Levels  of  significance  for  individual  comparisons  between  washout 
filter  types  are  shown  in  table  I. 

Figures  21  through  23  show  the  average  frequency  distributions  of  control  move- 
ments for  first  order  filters,  nonlinear  adaptive  filters  and  high  order  filters  with  tilt 
respectively. 

Comparisons  were  made  between  the  different  frequency  distrubitions  by  subtract- 
ing one  spectrum  from  another  and  limiting  the  difference  spectrum  to  zero  as  a minimum. 

The  difference  spectrum s for  cyclic  x control  movements  are  shown  in  the  follow- 


ing figures: 

First  Order  Washout  - Nonlinear  Adaptive  Washout  Figure  24 

Nonlinear  Adaptive  Washout  - First  Order  Washout  Figure  25 

First  Order  Washout  - High  Order  Washout  with  Tilt  Figure  26 

High  Order  Washout  with  Tilt  - First  Order  Washout  Figure  27 

Nonlinear  Adaptive  Washout  - High  Order  Washout  with  Tilt  Figure  28 

High  Order  Washout  with  Tilt  - Nonlinear  Adaptive  Washout  Figure  29 
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Control  Variance  For  Washout  Types  (Phase  III) 
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FIRST  ORDER  HIGH  ORDER  NONLINEAR 

WASHOUT  WASHOUT  ADAPTIVE 

WITH  TILT  WASHOUT 

FIGURE  18  - Body  Rate  Error  Variance  For  Washout  Types  (Phase  III) 
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TABLE  I 

LEVELS  OF  SIGNIFICANCE  FOR  PHASE  IH  OF  THE  EXPERIMENT 


Moasurc  of  Lt-vcl  of 

Effoctivcncss  Comparison  Significimco 

Control  Variance  First  Order  versus  High  Order  with  Tilt  . 13 

Control  Variance  High  Order  with  Tilt  versus  Adaptive  Nonlinear  .07 

Control  Variance  First  Order  versus  Adaptive  Nonlinear  .03 

Body  Rate  Error  First  Order  versus  High  Oi'der  with  Tilt  .08 

Variance 

Body  Rate  Error  High  Order  with  Tilt  versus  Adaptive  Nonlinear  . 21 

Variance- 

Body  Rate  Error  First  Order  versus  Adaptive  Nonlinear  . 03 

Variance 

Attitude  Error  First  Order  versus  High  Order  with  Tilt  . 07 

Variance- 

Attitude  Error  High  Order  with  Tilt  versus  Nonlinear  .Adaptive  . 10 

Variance- 

Attitude  Error  First  Order  versus  Nonlinear  Adaptive  . 04 

Variance- 

Position  Error  Comparison  of  all  three  . 75 

Variance- 


For  cyclic  y the  difference  spectrums  appear  in  the  following  figures; 

Figure  30 
Figure  31 
Figure  32 
Figeire  33 
Figure  34 
Figure  35 

For  collective  movements  the  difference  spectrums  appear  in  the  following  figures; 


First  Order  Washout  - Nonlinear  Adaptive  Washout  Figure  36 

Nonlinear  Adaptive  Washout  - First  Order  Washout  Figure  37 

First  Order  Washout  - High  Order  Washout  with  Tilt  Figure  38 

High  Order  Washout  with  Tilt  - First  Order  Washout  Figure  39 

Nonlinear  Adaptive  Washout  - High  Order  Washout  with  Tilt  Figure  40 

High  Order  Washout  with  Tilt  - Nonlinear  Ad.aptive  Washout  Figure  41 


First  Order  Washout  - Nonlinear  Adaptive  Washout 
Nonlinear  Adaptive  Washout  - First  Order  Washout 
First  Order  Washout  - High  Order  Washout  with  Tilt 
High  Order  Washout  with  Tilt  - First  Order  Washout 
Nonlinear  Adaptive  Washout  - High  Order  Washout  with  Tilt 
High  Order  Washout  with  Tilt  - Nonlinear  Adaptive  Washout 
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FIGUllE  21  - Frequency  Content  Of  Control  Movements  For  First  Order  Washout 
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FIGURE  21  - Frequency  Content  Of  Control  Movements  For  First  Order  Washout  (Con't) 
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FIGlMtf;  21  - Frequency  Content  Of  Control  Movements  For  First  Order  Washout  (Con't) 
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FIGURE  22  - Frequency  Content  Of  Control  Movements  For  Nonlinear  Adaptive  Washout 
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FIGURE  22  - Frequency  Content  Of  Control  Movements  For  Nonlinear  Adaptive  Washout  (Con't) 
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FIGURE  22  - Frequency  Content  Of  Control  Movements  For  Nonlinear  Adaptive  Washout  (Con't) 
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FIGURE  23  - Frequency  Content  Of  Control  Movements  For  High  Order  Washout  With  Tilt 
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FIGURE  23  - Frequency  Content  Of  Control  Movements  For  High  Order  Washout  With  Tilt  (Con’t) 
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FIGURE  24  - Cyclic  X Difference  Spectrum  For  First  Order  Washout  Minus  Nonlinear  Adaptive  Washout 
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FIGURE  25  - Cyclic  X Difference  Spectrum  For  Nonlinear  Adaptive  Washout  Minus  First  (mder  Washout 
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FIGURE  26  - Cyclic  X Difference  Spectrum  For  First  Order  Washout  Minus  High  Order  Washout  With  Tilt 
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FIGURE  27  - Cyclic  X Difference  Spectrum  For  High  Order  Washout  With  Tilt  Minus  First  Order  Washout 
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FIGURE  28  - Cyclic  X Difference  Spectrum  For  Nonlinear  Adaptive  Washout 
Minus  High  Order  Washout  With  Tilt 
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FIGURE  29  - Cyclic  X Difference  Spectrum  For  High  Order  Washout  With 
Tilt  Minus  Nonlinear  Adaptive  Washout 
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Minus  Nonlinear  Adaptive  Washout 
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FIGURE  31  - Cyclic  Y Difference  ^ectrum  For  Nonlinear  Adaptive  Washout 

Minus  First  Order  Washout 
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FIGURE  32  - Cyclic  Y Difference  Spectrum  For  First  Order  Washout 
Minus  High  Order  Washout  With  Tilt 
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FIGURE  33  - Cyclic  Y Difference  Spectrum  For  High  Order  Washout 
With  Tilt  Minus  First  Order  Washout 
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FIGURE  34  - Cyclic  Y Difference  Spectrum  For  Nonlinear  Adaptive  Washout 
Minus  High  Order  Washout  With  Tilt 
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Minus  Nonlinear  Adaptive  Washout 
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Minus  Nonlinear  Adaptive  Washout 
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FIGURE  37  - Collective  Difference  Spectrum  For  Nonlinear  Adaptive  Washout 

Minus  First  Order  Washout 
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FIGURE  38  - Collective  Difference  ^eetruni  For  First  Order  Washout 
Minus  High  Order  Washout  With  Tilt 
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/•7GURE  39  - Collective  Difference  Spectrum  For  High  Oi'der  Washout 
With  Tilt  Minus  First  Order  Washout 
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FIGURE  40  - Collective  Difference  Spectrum  For  Nonlinear  Adaptive  Washout 
Minus  High  Order  Washout  With  Tilt 
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Minus  Nonlinear  Adaptive  Washout 
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Figures  21,  22,  24,  and  25  indicate  that  the  frequency  content  for  cyclic  x move- 
ments are  nearly  the  same  for  first  order  filters  and  nonlinear  adaptive  filters.  There 
is  a tendency,  however,  for  the  energy  to  be  more  concentrated  in  the  lower  end  of  the 
spectrum  when  using  nonlinear  adaptive  filters  than  with  the  first  order  filters.  This 
is  an  indication  of  a tighter  control  loop  for  first  order  filters.  Apparently  the  pilots 
are  receiving  better  motion  cues  and  reacting  quicker  with  the  first  order  filters  than 
with  the  nonlinear  adaptive  filters.  The  same  trend  is  evident  for  both  cyclic  y and 
collective  control  movements  as  indicated  by  figures  30,  32,  36  and  37. 

Figures  21,  23,  26,  27,  32,  33,  38  and  29  indicate  that  more  energy  is  present  in 
the  spectrum  of  control  movements  when  first  order  filters  are  used  as  opposed  to  high 
order  filters  with  tilt.  It  is  also  evident  from  figures  26,  27,  and  33  that  enci'gy  is  con- 
centrated more  in  the  lower  end  of  the  spectrum  for  high  order  filters  with  tilt  than  for 
first  order  filters.  Both  results  indicate  a tighter  control  loop  when  first  order  filters 
are  used.  Figures  22,  23,  28,  29,  34,  35  and  41  are  comparisons  of  control  movement 
spectrums  for  nonlinear  adaptive  filters  and  high  order  filters  with  tilt.  There  seems 
to  be  a slightly  higher  energy  level  in  the  spectrum  for  cyclic  x when  using  nonlinear 
adaptive  filters  and  perhaps  for  cyclic  y also.  This  trend  is  reversed,  however,  for 
collective  movements.  Therefore  there  is  no  consistent  trend  for  differences  in  spec- 
trums when  using  the  nonlinear  adaptive  filters  versus  high  order  filters  with  tilt. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


From  the  pilot  comments,  results  of  the  statistical  tests  and  the  time  histories  of 
the  motion  algorithms,  it  is  apparent  that  the  axis  system  of  application  for  washout 
filters  has  no  effect  on  the  accuracy  with  which  pilots  fly  the  aircraft.  Indeed,  the  axis 
system  transformations  may  be  linearized  by  applying  the  small  angle  approximation 
without  jeopardizing  the  quality  of  the  motion. 

The  results  of  the  second  phase  also  indicate  that  high  order  washout  filters  without 
tilt  should  be  avoided  for  three-degree-of-freedom  motion  bases.  The  anomalous 
specific  force  cues  that  result  from  angle  washout  (a  characteristic  of  high  order  filters 
without  tilt)  are  probably  the  cause  for  pilot  induced  oscillations  and  resulting  higher 
error  variances  that  occur  when  this  t3qje  of  washout  is  used. 

The  fact  that  no  significant  difference  in  performance  can  be  observed  for  mns 
with  motion  as  opposed  to  nms  without  motion  is  probably  a result  of  the  flying  task 
used  in  this  e:q)eriment.  Visual  flight  tasks  are  better  for  measuring  the  relative  merit 
of  washout  circuits  than  instrument  flight  tasks.  Since  no  out-the-window  visuals  were 
available  for  this  simulation,  the  steering  commands  were  given  by  an  instrument  posi- 
tioned above  the  instrument  panel  in  an  attempt  to  require  the  pilot  to  set  up  a sc.ui  pat- 
tern similar  to  the  one  he  would  use  when  flying  a visual  task.  This  apparently  was  not 
totally  successfully,  because  thcpilots  Indicated  that  they  still  flew  as  if  it  were  ;in  in- 
strument flying  task.  Motion  is  not  as  important  a cue  when  flying  instruments  as  it  is 
when  flying  a visual  flight  task  such  as  carrier  landing  or  terrain  following.  For  the 
flight  task  used  in  this  experiment,  the  pilots  felt  that  motion  was  valuable  more  from 
a standpoint  of  adding  realism  to  the  simulation  rather  than  from  giving  them  additional 
cues  for  flying  the  aircraft. 

The  results  of  Phase  HI  show  conclusively  that  first  order  washout  filters  provided 
the  best  motion  cues  to  the  pilot  for  the  flying  task,  aircraft  type  and  motion  base  hard- 
ware used  in  this  experiment. 

Smaller  variances  for  control  movements,  body  rate  errors  and  attitude  errors 
can  be  attributed  to  a tightening  of  the  control  loop  when  using  first  order  filters  as 
evidenced  by  the  frequency  distributions  of  the  control  movements. 

Results  concerning  nonlinear  adaptive  washouts  and  high  order  washouts  with  tilt 
arc  inconclusive.  The  time  histories  in  Appendix  C indicate  that  first  order  filters  pro- 
vide better  motion  recovery  than  the  two  linear  filters  for  all  except  the  veitical  degree 
of  freodom.  However,  the  error  variances  indicate  that  nonlinear  adaptive  filters  did 
not  provide  the  best  motion  cues  in  the  e;q)erlment.  Reference^  offers  a possible 

3 

Graham,  D.  K. , Hickey,  L.  K.,  and  Vaden,  V.  C. ; "A  Rational  for  Moving  Base 
Flight  Simulation  and  a Preliminary  Statement  of  the  Motion  Requirements";  The 
Boeing  Company;  Oct  1968. 
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expl;uiation  of  this  seeming  contradiction.  Flight  in  turbulent  air  is  listed  as  an  ex- 
ample of  a type  of  motion  which  actually  hinders  the  pilot  when  flying  his  aircraft. 
Turbulence  was  passed  most  effectively  by  the  nonlinear  adaptive  filters,  less  effec- 
tively by  the  high  order  filters  with  tilt  and  least  effectively  by  the  first  order  filters. 

It  is  possible  that  confusion  caused  by  the  turbulence  was  responsible  for  the  poor 
performance  of  nonlinear  adaptive  washouts  in  this  experiment. 

Addition  of  tilt  to  the  high  order  washouts  did  not  appreciably  improve  the  specific 
force  recovery  and  is  therefore  not  recommended  for  simulation  of  helicopter  motion 
on  three  degree-of- freedom  motion  bases.  The  frequency  of  specific  forces  encoim- 
tered  in  helicopter  motion  is  too  high  to  be  passed  by  the  tilt  shaping  filter.  There- 
fore, the  anomalous  rotation  cues  caused  by  the  high  order  rotational  washout  c;innot 
be  counteracted  by  tilt. 

Once  again  the  pilots  indicated  that  nonlinear  adaptive  filters  produced  the  most 
realistic  motion  because  the  turbulence  was  most  apparent  when  using  that  kind  of 
filter.  However,  they  felt  that  all  three  filter  types  (first  order,  nonlinear  adaptive 
and  high  order  with  tilt)  produced  good  quality  motion  simulation. 

Position  error  variance  was  not  found  to  be  a sensitive  measure  of  effectiveness 
for  comparing  motion  algorithms.  It  is  suggested  that  future  studies  of  this  sort  could 
benefit  by  collecting  data  on  body  rates  and  accelerations,  linear  velocities  and  accel- 
erations, and  control  movements. 
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APPENDIX  A 

DIGITAL  TURBULENCE  MODEL 


The  vertical  con\ponent  of  turbulent  air  as  seen  by  a moving  aircraft  has  a pov  er 
spectral  density  given  by: 


<1.^  (co)  = Wg2 


1^ 

rr 


2 


2 


where: 

2 

W is  the  mean  square  value  of  the  vertical  gust  - a simple  measure  of 
^ the  intensity  of  the  turbulence. 

L is  the  scale  of  the  turbulence  - a measure  of  the  average  eddy  size,  it 
indicates  the  frequency  at  which  the  spectrum  breaks  from  a flat  spectrum 
to  one  that  decreases  rapidly. 

V is  the  forward  velocity  of  the  aircraft. 

A 

Reference'll  a report  dealing  with  the  simulation  of  atmospheric  turbulence,  indi- 
cates that  vertical  gust  turbulence  may  be  simulated  by  passing  white  Gaussian  noise 
through  a filter  having  the  following  transfer  function: 


Experimental  data  indicates  the  L is  of  the  order  of  1000  feet.  With  respect  to  the 
intensity  of  the  turbulence,  a root  mean  square  value  of  15  feet/sec  was  used  in  this 
simulation.  The  Intensity  of  turbulence  is  known  to  be  subject  to  random  fluctuations. 


<1 

l{.  lU  rkowitz  and  J.  Loftus;  "Simulation  of  Atmospheric  Turbulence";  U.  S.  Naval 
Air  Development  Center,  Report  ACLTN  No.  62-147-02;  Feb  1962. 
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At  low  altitude,  root  mean  square  values  of  vertical  gust  velocities  seem  to  have  a 
Rayleigh  probability  distribution.  No  data  is  available  concerning  the  rate  of  change 
of  the  intensity  of  turbulence,  however,  an  average  period  of  30  seconds  seems  reason- 
able and  was  used  in  this  simulation. 


Therefore,  the  root  mean  square  vertical  gust  velocity  was  generated  by  first 
sampling  a Itayleigh  distribution  with  a mean  of  15  ft/sec  and  passing  the  random  vari- 
abb  through  a low  pass  filter  with  the  following  characteristics: 

1 

30  S ^ 1 

2 

The  value  thus  derived  is  Wg  and  is  used  to  sample  a Gaussian  distribution.  The 
random  variable  derived  from  this  sample  is  then  passed  through  the  filter  mentioned 
previously  which  reproduces  the  spectrum  of  atmospheric  vertical  gust  turbulence. 

iteference^,  in  the  chapter  on  flight  in  turbulent  air,  suggests  that  the  veitical  gust 
velocity  be  input  to  the  aircraft  model  in  the  following  manner.  Three  aircraft  state 
variables  are  most  affected  by  vertical  gusts:  the  angle  of  attack,  roll  rate  and  pitch 
rate.  Perturbations  to  these  quantities  caused  by  vertical  gust  can  be  generated  as 
follows  and  added  to  the  aircraft  simulation: 


a 


K 


P 


-W  /V 

g 


A 


-W  /V 
g 

W /V 
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A 

A 


where: 


a 

g 


is  the  angle  of  attack  perturbation  due  to  vertical  gusts. 


P 

g 


is  the  roll  rate  perturbation  due  to  vertical  gusts. 


is  the  pitch  rate  perturbation  due  to  vertical  gusts. 


W 

g 


is  the  Instantaneous  vertical  gust  velocity. 


B.  Etkin;  "Dynamics  of  Flight  — Stability  and  Control";  John  Wiley  and  Sons,  Inc. ; 
1950. 


A-2 


APPENDIX  B 

MOTION  BASE  HYDRAULIC  MAINTENANCE 


NADC-77306-20 


APPENDIX  B 

MOTION  BASE  HYDRAULIC  MAINTENANCE 


The  motion  base  drive  investigation  required  the  integration  of  generic  helicopter 
flight  motion  equations  with  the  Link  three-degree-of-freedom  motion  base  system. 

The  following  description  of  motion  base  hydraulic  maintenance  is  an  excei'pt  from  a 
technical  memorandum  written  by  Mr.  T.  Thomas  of  the  Performance  Simulation 
Division. 

The  motion  system  was  found  to  require  some  hydraulic  refurbishment  in  t)rder 
that  the  system  could  be  fully  utilized.  The  following  S3rmptoms  were  noted: 

1.  The  hydraulic  pressure  on  the  sequence  valve  was  down  from  3000  psi  to  about 
1400  psi. 

2.  The  motion  of  the  three  actuator  legs  was  not  fully  responsive.  If  an  input  sig- 
nal was  given  to  an  actuator  leg,  an  unreasonable  delay  time  was  noted  before  the  actu- 
ator leg  was  driven  to  the  commanded  position.  For  example  the  right  actuator  leg 
drive  at  a maximum  speed  of  about  3.5  inches  per  second. 

3.  The  motion  of  the  three  actuator  legs  did  not  occur  at  the  same  rate. 

4.  The  motion  of  the  three  actuator  legs  was  jerky. 

The  Fluid  System  Section  (code  30211)  of  the  Air  Vehicle  Technology  Department 
was  consulted  regarding  these  hydraulic  problems.  Mr.  J.  F.  Ohlson  (Ext  2834)  heads 
this  section.  Their  technicians  inspected  the  motion  system  and  made  several  sug- 
gestions regarding  the  motion  systems  condition  and  refurbishment.  The  main  relief 
valve  (tag  number  10  on  the  unit)  was  warm  during  the  motion  base  operation  which 
meant  that  hydraulic  fluid  was  being  spilled  out  of  this  relief  valve  and  not  through  the 
required  sequence  valve  (tag  number  7 on  the  unit).  The  main  relief  valve  was  adjusted 
but  the  system  pressure  could  not  be  raised  above  1500  psi  utilizing  just  this  relief  valve 
adjustment.  It  was  thought  that  this  relief  valve  might  be  remaining  open  and  therefore 
required  cleaning  and  the  replacement  of  the  valves  O rings. 

Mr.  Colson  agreed  to  allow  one  of  his  hydraulic  technicians  to  spend  one  week 
analyzing  the  motion  systems  problems  and  refurbishing  the  system  where  necessary. 
The  technician  assigned  to  the  job  was  Mr.  M.  Brandt  (Ext.  2062). 

The  following  things  were  done  to  the  motion  system  In  an  attempt  to  improve  its 
peiTormance. 

1.  The  main  relief  valve  (No.  10)  was  taken  apart,  cleaned,  inspected  and  the 
O rings  replaced. 
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2.  The  right,  left  and  rear  cylinder  filters  were  replaced.  It  was  observed  that 
the  rear  cylinder  bowl  was  slightly  dirty  and  the  old  filters  were  discolored.  The  filter 
number  and  spec  is  9C4330-00-277-3274,  QPL  elements,  10  micron  25  abs,  ANG235-4A. 

3.  The  hydraulic  fluid  reservoir  holding  about  50  gallons  of  oil  was  pumped  out. 

The  oil  utilized  by  the  Singer  Company  was  identified  as  micronic  745  oil.  The  military 
uses  MIL-H-5606C  hydraulic  fluid  (Stock  No.  9G-9150-00-223-4134)  in  their  aircraft. 
Sixty-five  gallons  of  this  5606C  t5rpe  oil  was  obtained  from  the  Center's  stock.  The 
5G0GC  and  the  745  type  oils  are  compatible  according  to  Mr.  Larry  Almy  of  the  Singer 
Company.  (His  telephone  number  is  1-607-772-3363).  The  oil  removed  from  the 
reservoir  smelled  burnt  which  might  indicate  that  the  water  cooling  system  might  not 

be  adequate. 

4.  A new  drain  spigot  was  placed  on  the  bottom  of  the  reservoir  to  facilitate  future 
removal  of  oil  from  the  reservoir. 

5.  The  front  plate  on  the  reservoir  was  removed  and  the  metal  strainer  taken  out 
of  the  reservoir.  The  strainer  was  cleaned  and  put  back  into  the  reservoir. 

6.  The  old  745  type  oil  was  flushed  out  of  the  hydraulic  lines  leading  from  the  pump 
assembly  to  the  motion  system  and  back  to  the  reservoir.  This  was  accomplished  by 
removing  the  main  return  line  from  the  motion  system,  turning  on  the  pump  and  motion 
systems  and  forcing  the  fluid  out  of  the  main  return  line  into  a 55  gallon  empty  drum. 
Before  turning  on  the pump  and  motion  system,  the  reservoir  was  filled  about  half  way 
with  new  5606C  type  fluid.  It  was  observed  that  one  inch  of  the  reservoir  indicator  tube 
is  equal  to  about  six  gallons  of  hydraulic  fluid. 

The  pump  and  motion  system  were  then  turned  on.  As  the  motion  base  moved  up- 
ward the  old  745  oil  flowed  into  the  55  gallon  waste  drum.  Since  the  745  oil  has  a 
brownish  yellow  color  and  the  new  5606C  oil  has  a redish  color,  it  was  easy  to  observe 
when  all  the  old  oil  was  flushed  out  of  the  motion  base  hydraulic  lines. 

Each  intake  line  to  the  right,  left  and  rear  actuators  was  cracked  open  separately 
to  insure  that  any  air  in  the  lines  would  be  removed.  With  the  motion  base  on,  it  was 
observed  that  in  addition  to  a slight  amount  of  hydraulic  fluid  leaking  out  of  the  intake 
lines,  bubbles  also  were  observed.  Before  the  air  was  removed  from  the  actuator 
intake  lines  the  right  actuator  would  not  move.  This  very  well  may  have  been  caused 
by  a vapor  lock  condition  in  the  right  actuators  solenoid  system.  More  hydraulic  fluid 
was  added  to  the  reservoir  until  the  glass  measuring  gage  was  full.  The  total  system 
requires  about  65  gallons  of  hydraulic  fluid. 

7.  The  hydraulic  filters  were  again  replaced  with  new  filters. 
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8.  The  hydraulic  valve  system  was  then  adjusted  to  see  if  the  system  pressure 
could  be  increased  to  3000  psi.  The  Singer  Company  Hydraulic  Troubleshooting  Manual 
was  utilized  for  this  operation.  The  steps  described  in  the  manual  were  followed  to 
increase  the  system  pressure  to  3000  psi.  The  procedure  in  the  manual  is  fairly  good 
but  it  helps  to  know  something  about  hydraulic  systems  to  supplement  some  of  the  manu- 
als steps.  For  example,  step  5 describes  the  method  of  increasing  the  system  pres- 
sure to  3400  psi  in  order  that  the  main  relief  valve  (No.  10)  can  be  adjusted  to  its  proper 
pressure.  This  is  accomplished  by  adjusting  the  sequence  valve  (No.  5).  Adjusting  the 
sequence  valve  (No.  5)  to  its  maximum  only  increased  the  system  pressure  to  2000  psi. 
Adjusting  the  main  relief  valve  (No.  10)  raised  the  system  pressure  to  3400  psi.  The 
sequence  valve  (No.  5)  was  then  adjusted  back  to  3000  psi  (Step  No.  8 in  the  manual). 

9.  A thermometer  system  was  added  to  the  reservoir  so  that  the  oil  temperature 
can  be  observed  during  continual  motion  system  operation.  This  will  indicate  the  ade- 
quacy of  the  motion  systems  cooling  system. 

10.  The  motion  system  required  a new  solenoid  valve  in  the  left  actuator  leg. 

A high  pitch  squeal  was  detected  whenever  the  left  actuator  was  driven.  Tne  solenoid 
valve  was  replaced  with  an  old  Moog  type  solenoid  valve  that  had  been  used  in  the  mo- 
tion system  a couple  of  years  ago.  Actually  the  right  solenoid  valve  was  put  on  the  left 
actuator  leg  and  the  old  Moog  solenoid  valve  was  put  on  the  right  actuator  leg. 

The  motion  system  was  checked  to  see  if  there  was  any  improvement  in  perform- 
ance of  the  motion  system  after  refurbishment.  The  system  showed  only  a slight  im- 
provement in  performance  and  the  jerky  motion  of  the  actuator  legs  was  still  detected. 

The  motion  systems  accumulator  pressure  is  required  to  be  set  at  1400  psi.  The 
valve  on  the  top  of  the  accumulator  was  found  to  have  a leak  and  consequently  the  ac- 
cumulator pressure  was  down  to  just  about  0 psi.  The  valve  was  replaced  and  the 
accumulator  was  recharged  with  nitrogen  to  the  required  1400  psi.  Mr.  Phil  Porter 
(Ext.  2638)  of  the  Equip/Inst  Shop  Branch  is  the  technician  that  charged  the  accumulator. 

The  motion  system  was  again  checked  to  see  if  there  was  any  improvement  in  per- 
formance. The  system  now  showed  a marked  improvement  in  performance  and  the  jerky 
motion  of  the  actuator  legs  was  eliminated.  The  right  and  rear  legs  were  driven  at  a 
maximum  velocity  of  about  14  inches  per  second.  There  was  still  some  slight  perform- 
ance problems  with  the  right  actuator  leg.  The  rear  actuatcr  solenoid  valve  was  put 
on  the  right  actuator  leg  and  a different  old  Moog  t}rpe  solenoid  valve  was  put  on  the  rear 
actuator  leg.  The  right  actuator  leg  appears  to  operate  fairly  well  but  the  rear  actu- 
ator leg  drove  with  a jerky  motion.  A third  old  solenoid  valve  was  put  on  the  rear 
actuator  leg  and  this  time  all  three  actuator  legs  appeared  to  have  improved  perform- 
ance. 
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rhf  motion  system  was  again  checked  to  see  if  there  was  any  improvement  in  per- 
formance. The  system  showed  acceptable  performance  for  all  three  actuator  legs. 

The  maximum  velocity  of  all  three  actuator  legs  was  about  14  inches  per  second. 

The  final  step  in  the  refurbishment  of  the  motion  base  was  the  adjustment  of  the 
electronic  gains  in  the  actuator  leg  servo  systems. 

Optimal  adjustment  of  the  electronic  gains  resulted  in  actuator  leg  responses 
similar  to  a second  order  system.  For  the  front  legs,  the  gains  could  be  adjusted  to 
achieve  a 3 dB  bandwidth  of  1 Hz.  The  rear  leg  experienced  high  frequency  oscilla- 
tions when  the  3 dB  bandwidth  was  increased  to  . 75  Hz  or  greater. 
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WASHOUT  CIRCUIT  DESCRIPTIONS 


Each  of  the  twelve  motion  algorithms  investigated  in  the  second  phase  of  the  ex- 
periment are  described  herein.  Included  in  each  description  is  a block  diagram  illus- 
trating the  computational  flow  for  the  algorithm. 

Euler  transformation  [e],  axis  transformation  [f]  and  and  simulator 

transformation  from  pitch,  roll  and  heave  command  to  leg  commands  are  detailed  at 
the  end  of  this  appendix. 

FIRST  ORDER  WASHOUT  (BODY  AXIS) 

Figure  C-1  is  a block  diagram  of  a first  order  body  axis  washout  circuit  for  a three 
degree-of-freedom  motion  base.  The  aircraft  body  axis  roll  and  pitch  rates  (Pj^andQ  :^) 
represent  the  desired  angular  motion  for  the  simulator.  The  P^  and  Qa  signals  are 
passed  through  first  order  washout  filters  with  gains  Kp  and  Kq  both  equal  to  0.5  and 
frequency  breakpoints  Wp  and  Wq  equal  to  0.5  rate/sec.  The  output  of  the  washout 
filters,  Pi  and  Qi,  represent  the  achievable  body  axis  roll  and  pitch  rates  for  the 
simulator. 

Note  that  Rj,  the  body  axis  yaw  rate,  must  be  calculated  from  and  the  timgent 
of  <l>  jvi,  the  Euler  roll  angle  of  the  simulator,  to  assure  that  ♦ ]yj,  the  Euler  yaw  angle 
of  the  simulator,  will  always  be  zero. 

Pi,  Qi  and  Ri,  the  achievable  body  axis  rates  for  the  simulator,  are  resolved  into 
Euler  angular  rates  based  on  the  current  simulator  attitude.  As  mentioned  above, 'I'"  ^ 
will  always  be  zero.  The  Euler  roll  and  pitch  rates,  0 and  are  integrated  to 
produce  the  roll  and  pitch  angle  commands,  <t>  and  0j^. 

The  Z axis  linear  drive  represents  the  third  degree  of  freedom.  Wp^  is  the  air- 
craft vertical  acceleration  and  represents  the  desired  vertical  acceleration  for  the 
simulator.  The  acceleration  due  to  gravity,  G,  is  added  to  Wp^  to  insure  that  only 
deviations  of  vertical  acceleration  from  one  G are  washed  out.  The  vertical  accelera- 
tion deviation  is  passed  through  a third  oi’der  washout  filter  with  constants  C32  = 1.  77, 
C4/  - 1.50  and  C52  = .125.  The  full  travel  of  the  simulator  Z axis  is  very  limited, 
thus  necessitating  the  use  of  a third  order  washout.  The  third  order  washout  filter  will 
wash  out  position  as  well  as  acceleration,  thereby  limiting  the  travel  of  the  simulator 
more  effectively  than  a lower  order  washout  filter.  Z2,  the  achievable  vortical  ac- 
celeration deviation  from  one  G,  is  then  doubly  integrated  and  multiplied  by  a gain 
factor  of  K2  = . 1 to  produce  Z2,  the  simulator  Z axis  command. 
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- First  Order  Body  Axis  Washout  Circuit 
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The  roll,  pitch  and  heave  commands  ® M the*'  transformed 

to  simulator  leg  commands.  The  equations  for  performing  this  transformation  are 
listed  at  the  end  of  this  appendix. 

FIRST  ORDER  WASHOUT  (INERTIAL  AXIS) 

Figure  C-2  shows  a block  diagram  of  a first  order  inertial  axis  washout  circuit. 
The  desired  body  rates,  and  are  first  transformed  to  desired  Euler  angular 
rates,  ^ ^ and  q j.  Once  again  R^  is  calculated  to  insure  U^at^’  i will  always  be 
zero.  The  desirable  Euler  angular  rates  are  passed  throu  first  order  washout  filters 
with  characteristics  Kp  = Kj^  = 0.5  and  Wp  = = 0. 5 to  obtain  the  achievable  rates 

and  ’©M.  Integration  of  the  achievable  rates  produces  the  simulator  roll  and 
pitch  command,  ^ and  6 The  vertical  acceleration  is  identical  to  that  of  the  first 
order  body  axis  washout  circuitry. 

FIRST  ORDER  WASHOUT  (LINEAR  TRANSFORMATION) 

Figure  C-3  is  a block  diagram  of  a linearized  first  order  washout  circuit.  The 
desired  roll  and  pitch  rates,  and  Q^,  are  passed  through  first  order  washout 
filters  cascaded  with  integrations  to  yield  the  simulator  roll  and  pitch  commands,  4 m 
and  6 The  vertical  acceleration  is  once  again  identical  to  that  for  the  other  first 
order  washout  circuits. 

THIRD  ORDER  WASHOUT  CIRCUIT  (BODY  AXIS) 

Figure  C-4  is  a block  diagram  of  a third  order  body  axis  washout  circuit.  Py^ 
and  Qy^  are  passed  through  third  order  washout  filters  with  characteristics:  Kp  = 

Kp  = . 6,  W2P  = ^2R  “ • ^3P  ~ ^3R  ® *3  and  Zgp  = = 1.0.  The  output  of  the 

filters.  Pi  and  Qi  are  transformed  to  Euler  angular  rates,  and  and  inte- 

grated to  yield  the  simulator  roll  and  pitch  commands,  ♦ m and©  j^.  The  vertical  ac- 
celeration for  third  order  washout  circuits  is  handled  in  the  manner  as  the  first  order 
washout  circuits. 

THIRD  ORDER  WASHOUT  CIRCUIT  (INERTIAL  AXIS) 

Figure  C-5  is  a block  diagram  of  a third  order  inertial  axis  washout  circuit. 

Py^  and  Qy^  are  first  transformed  to  the  Euler  angular  rates,  *(^i  and  ’©j^,  before 
being  passed  through  third  order  washout  filters  cascaded  with  integrations  to  yield 
the  simulator  roll  and  pitch  commands. 

THIRD  ORDER  WASHOUT  CIRCUIT  (LINEAR  TRANSFORMATION) 

A linearized  third  order  washout  circuit  is  shown  in  figure  C-6.  P^  and  Qy^  are 
passed  directly  through  third  order  washout  filters  cascaded  with  integrations  to  yield 
the  simulator  roll  and  pitch  commands. 
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FIGURE  C-2  - First  Order  Inertial  Axis  Washout  Circuit 
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- Third  Order  Body  Axis  Washout  Circuit 
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[•’IGURE  C-5  - Third  Order  Inertial  Axis  Washout  Circuit 


C-6  - Linearized  Third  Order  Washout  Circuit 
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HIGH  ORDER  WASHOUT  CIRCUIT  WITH  TILT  (BODY  AXIS) 

Figure  C-7  is  a block  diagram  of  a high  order  body  axis  washout  circuit  with  tilt. 

and  are  passed  through  first  order  washout  filters  with  characteristics  Kp  = 

Kq  = . 3028  and  Wp  = Wq  = .4  to  produce  and  the  contribution  from  the  aircraft 
rotational  degrees  of  freedom.  These  rotation  rates  are  summed  with  P2  and  Q2,  the 
lateral  and  longitudinal  tilt  rates  to  produce  P and  Q,  the  total  body  axis  roll  and  pitch 
rates  for  the  simulator.  R is  calculated  to  insure  a of  zero  and  P,  Q and  R are 
transformed  to  0 m "®M*  integrations  then  yield  simulator  roll  and  pitch 

commands,  0 and 

• • 

The  aircraft  body  axis  longitudinal,  lateral  and  vertical  accelerations  (Up^,  Vp^ 
and  Wpp)  represent  the  desired  accelerations  and  are  used  to  generate  the  tilt  rates. 

UpQ  and  Vp^  are  limited  to  ±.  IG's  before  processing  to  insure  that  simulator  tilt  angles 
remain  at  reasonable  levels.  The  vertical  acceleration  for  the  high  order  tilt  algorithms 
is  the  same  as  for  the  first  and  third  order  washout  circuits.  After  the  deviations  from 
one  G have  been  washed  out,  -G  is  added  back  in  and  the  three  modified  accelerations, 
limited  Upo  and  Vp^  and  washed  out  WpQ,  are  transformed  to  inertial  axes  based  on  the 
current  simulator  attitude.  The  inertial  Z component,  Zjyjo*  is  discarded,  and  the 
longitudinal  and  lateral  components,  Xj^q  and  represent  the  difference  between 

the  actual  simulator  accelerations  and  the  desired  simulator  accelerations.  The  cross 
product  and  transformation  back  to  body  axes  yields  the  desired  longitudinal  and  lateral 
tilt  rates,  Pp  and  Qr,  and  Rx  is  discarded.  Px  and  Q-j.  are  passed  through  a shaping 
network  designed  to  minimize  rotational  anomalies  and  the  outputs,  body  axis  tilt  rates 
P2  and  Q2,  are  summed  with  the  rotational  channel  rates.  Pi  and  Q2.  Characteristics 
of  the  shaping  filter  are  as  follows: 


W = W 
IP  IR 


.004 


^8P  = ^8R 


.707 


= W„„ 

3P  3R 


.2 


The  tilt  rates  rotate  the  cockpit  to  null  out  Xjviq  and  Yj^iq-  When  these  two  inertial 
components  are  zero,  th^  the  simulator  longitudinal  and  lateral  body  axis  accelerations 


are  equal  to  the  limited  Up^  and  Vpo. 
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FIGURE  C-7  - High  Order  Body  Axis  Washout  Circuit  With  Tilt 
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HIGH  ORDER  WASHOUT  CmCUIT  WITH  TILT  (INERTIAL  AXIS) 

A block  diagram  of  a high  order  inertial  axes  washout  circuit  with  tilt  is  shown  in 
figure  C-8.  P^,  and  R^  (calculated  to  insure  ’'k  = 0),  are  transferred  to  Euler 

angular  rates  and  passed  through  first  order  filters  to  yield  2^  and  0^^.  These  rates 
represent  the  achievable  rotational  motion  derived  from  the  aircraft  rotational  degrees 
of  freedom.  The  Euler  tilt  rates,  i 2 ^nd  *62,  are  summed  with  ^ 1 and  to  form 
]yj  and  *0  jyj  the  total  rotational  motion  of  the  simulator.  ^ jyj  and  0^^  are  then  inte- 
grated to  yield  the  simulator  roll  and  pitch  command,  jyj  and  G^. 

For  the  tilt  channel,  the  error  signals,  XjyjQ  and  are  derived  as  for  the  body 

axis  tilt  algorithm:  limited  Up^  and  VpQ  and  washed  out  Wp^  are  transformed  to 
inertial  axes.  A desired  tilt  rate  vector  resolved  in  inertial  axes  is  generated  by 
taking  the  cross  product  of  the  vector 


X:' 

0 

J 

0 

^Mo 

and 

0 

s 

0 

1  

-1 

The  desired  x and  y body  axis  tilt  rates,  Pq^.  and  Qq.,  are  obtained  by  transforming  the 
desired  inertial  tilt  rates  to  body  axes.  Pq.  and  Qq,  are  passed  through  shaping  filters 
to  produce  P2  and  Q2,  the  body  axis  tilt  rates,  and  then  transformed  to  the  Euler  tilt 
rates,  0 2 '®2'  summed  with  the  rates  from  the  rotational  channel. 

HIGH  ORDER  WASHOUT  CIRCUIT  WITH  TILT  (LINEAR  TRANSFORMATIONS) 

A block  diagram  of  a linearized  high  order  washout  circuit  with  tilt  is  shown  in 
figure  C-9.  It  is  similar  to  the  other  tilt  algorithms  except  that  rotational  transforma- 
tions are  ignored  and  the  transformation  of  body  axis  accelerations  to  inertial  axes  is 
linearized  using  the  small  angle  approximation. 

NONLINEAR  ADAPTIVE  WASHOUT  CIRCUIT  (BODY  AXIS) 

The  block  diagram  in  figure  C-10  shows  the  computational  flow  for  a non-linear 
adaptive  body  axis  washout  circuit.  Adaptive  parameters  P^,  P0and  Pg  are  gen- 
erated to  minimize  the  cost  functions: 

= 1/2  (P^-P)^  + By/2  ( ypdt)^  + Cy/2(P)^  ^ Wy  (Y^^)^ 

J0=  1/2  (Q^-Q)^  + B^/2  ijQdt)^  + C^/2(Q)^  + 

Jz  = 1/2  (W^  + g -i/  + + D^/2 
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FIGURE  C-8  - High  Order  Inertial  Axis  Washout  Circuit  With  Tilt 


■.0-A053  630 
UNCLASSIFIED 


NAVAL  AIR  development  CENTER  WARMINSTER  PA  SYSTEMS  O—ETC  F/6  5/9 
investigation  of  MOTION  B»SE  DRIVE  TECHNIQUES* (U) 

MAR  76  M E BITNER 

NAOC-77306-a0  NL 


20*2 

*ftisa830 

B 

< i 

1 ) ; 

^ ► • 

1 * 

\ 

j * ^ 

’ n , 

«',r" 

• 1 

» *, 

i 

i-  ^ 

■,  <; 

> 

1 ' ''  ' 

i 

i V . - , 

• • • > 

[fit 

t 1 » » 

1 ) 

m 

1 1 • • 

m 

> M »’ 

V 

1 1 
) ) 

:ri!' 

• • • • 

M » » 

!! 

i » 

; f ; 
a.' 

t $ * 9 

ut 

» 1 I 1 

i 1 

1 \ 

■Vi 

» ) 1 > 

& 

' ''  '''  ; 

■ ■ • ■ 

s 

» ) 

m 
• 1 • 1 

► i * t ' 

SI 

) ) 

lit! 

■ s 

t » } » 

a 

> » 

• 1 • ■ 

lii"' 

( h t » 

i 

■1*' 

• > 1 < 

It 

i i t i 

ii 

) » 

!i 

■ ■ \ 

• 1 • 1 

M 

M > r 

n 

1 » 

• • • • 1 

> ) 1 > 

I 

» » 

IT* 

1 ( 

1 1 t 

END 

DATE 

FIIMEO 

6 -78 

DOC 

I 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS  l%J  A 


NADC-7  7306-20 


C-13 


NADC-7730(i-20 


C-14 


FIGURE  C-10  - Nonlinear  Adaptive  Body  Axis  Washout  Circuit 
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FIGURE  C-10  - Nonlinear  Adaptive  Body  Axis  Washout  Circuit  (Con't) 


NADC-77306-20 


C-IO  - Nonlinear  Adaptive  Body  Axis  Washout  Circuit  (Con't) 
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Equations  for  P^,  P0  and  are  derived  by  applying  the  method  of  continuous  steepest 
descent; 


J. 

P.  = -K.  V- 

1 1 P. 

9 ^ 


For  and  J2  this  becomes: 


aY 


-W„  Y 


M 


Y ‘M  gP 


'■z  ' ’^z 


I (W^+g-Z^)^  -c„  Z„  ^ -D„  Z ^ ^ 


2'  aPz 


z 2 3P^ 


z 2 aP 


Z 


] 


The  adaptive  filters  are  defined  as  follows: 

P = (P^  Pa  -Ey/ 

^ “®xy ^ 

Z2  = p^  (W^  +g)  -A^  Z2  -B^  Z^. 


The  state  sensitivity  equations  are  then  derived  from  the  definition  of  the  adaptive 
filters  and  the  body  axis  to  inertial  axis  transformation: 


5Z 

‘"a  -*z  >r: 


-B 


'I2 

Z dP. 
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d 

a 


e 


„ SIN  * +v„  SIN  „COS 

^“limited  ^“limited  ® 


-w^  cos*,  cose  \x  •-§-  dt 

^LIMITED  -^e 


d P* 


SIN*,-W„^  COS*  ,JX/ dt 

^°LIMITED  ^ ^LIMITED  ^/  •'  ^ 


The  output  of  the  rotational  adaptive  filters,  P and  Q,  are  transformed  to  Euler 
angular  rates  and  integrated  to  form  the  simulator  roil  and  pitch  commands,  0 and 

Bm- 


The  output  of  the  vertical  chaimel  adaptive  filter  is  the  simulator  heave  command. 


Z2- 


Washout  parameter  settings  are: 


Bx  = By  = C.25 


W =W  = .0001 

^ X 


= .09 

LIMITED  to  -29  and  -35  ft/sec^ 
Po 

U„  and  LIMITED  to  ±5  ft/sec^ 
Po  Po 
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Initial  conditions  for  the  state  variables  are  all  zero  except: 


NONLINEAR  ADAPTIVE  WASHOUT  CIRCUIT  (INERTIAL  AXIS) 

Shown  in  figure  C-11  is  a block  diagram  of  an  inertial  axis  washout  circuit.  Deriva- 
tion of  the  equations  is  analogous  to  that  for  the  body  axis  adaptive  washout  equations . 

For  the  inertial  axis  washout,  the  adaptive  filters  are  applied  to  Euler  angular  rates 
rather  than  body  axis  angular  rates.  The  vertical  channel  is  identical  to  that  for  the 
body  axis  adaptive  washout. 


NONLINEAR  ADAPTIVE  WASHOUT  CIRCUIT  (LINEAR  TRANSFORMATION) 

Figure  C-12  shows  a block  diagram  for  a nonlinear  adaptive  washout  circuit  for 
which  the  transformations  have  been  linearized.  The  vertical  channel  is  once  again 
the  same  as  for  the  body  and  inertial  axis  washout  circuits. 


Figures  C-13  through  C-24  show  the  response  of  each  of  the  washout  circuits  to 
the  same  maneuver.  The  maneuver  was  a turn  to  the  left  followed  by  a climb. 


The  Euler  transformation  matrix  for  a link  three-degree-of-freedom  motion  base 
is  as  follows: 


[E]  = 


0 


SIN0 

cos* 

TANe 

TANG 

COS^ 

-SING 

SIN* 

CCS* 

cose 

COSG 

where 


* 

p 

e 

= [e] 

Q 

R 
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FIGURE  C-11  - Nonlinear  Adaptive  Inertial  Axis  Washout  Circuit 
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FIGURE  C-11  - Nonlinear  Adaptive  Inertial  Axis  Washout  Circuit  (Con't) 
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FIGURE  C-11  - Nonlinear  Adaptive  Inertial  Axis  Washout  Circuit  (Con't^ 


►©Qa-iex  ®iii> 


NADC-77306-20 


a 

o 

u 

u 


s 

o 


0) 

> 

I 

<; 

0) 

.9 

0 

'z, 

01 
rH 

u 

H 

a 


C-24 


With  Linear  Transformations  (Con't) 
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FIGURE  C-12  - Nonlinear  Adaptive  Washout  Circuit  With  Linear  Transformations  (Con't) 
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FIGURE  C-13  - First  Order  Body  Axis  Washout  Circuit  Time  History 
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FIGURE  C-13  - First  Order  Body  Axis  Washout  Circuit 
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FIGURE  C-14  - First  Order  Inertial  Axis  Washout  Circuit  Time  History 
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FIGURE  C-14  - First  Order  Inertial  Axis  Washout  Circuit  Time  History  (Con't) 


FIGURE  C-14  - First  Order  Inertial  Axis  Washout  Circuit  Time  History  (Con't) 
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FIGURE  C-15  - Linearized  First  Order  Washout  Circuit  Time  History 
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FIGURE  C- 
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FIGURE  C-16  - Third  Order  Body  Axis  Washout  Circuit  Time  Ifistoiy 
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FIGURE  C-16  - Third  Order  Body  Axis  Washout  Circuit  Time  Histoir 
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FIGURE  C-IC  - Third  Ordor  Body  .\xis  Washout  Circuit  Time  History  (Con't) 
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FIGURE  C-17  - Third  Order  Inertial  Axis  Washout  Circuit  Time  History 
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FIGURE  C-17  - Third  Order  Inertial  Axis  Washout  Circuit  Time  History  (Con't) 
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FIGURE  C-18  - Linearized  Third  Order  Washout  Circuit  Time  History 
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FIGURE 
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FIC?URE  C-19  - Time  History  For  High  Order  Body  Axis  Washout  Circuit  With  Tilt 
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FIGURE  C-19  - Time  History 
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FIGURE  C-19  - Time  History  For  High  Order  Body  Axis  Washout  Circuit 
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FIGURE  C-20  - Time  History  For  High  Order  Inertial  ^\xis  Washout  Circuit  With  Tilt 


FIGURE  C-20  - Time  History  For  High  Order  Inertial  Axis  Washout  Circuit  With 
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FIGURE  C-20  - Time  History  For  High  Order  Inertial  Axis  Washout  Cireuit  With  Tilt  (Con't) 


FIGURE  C-21  - Time  History  For  Linearized  High  Order  Washout  Circuit  With  Tilt 
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FIGURE  C-21  - Time  History  For  Linearized  High  Orde: 
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FIGURE  C-22  - Time  History  For  Nonlinear  Adaptive  Rody  Axis  Washout  Circuit 


FIGURE  C-22  - Time  History  For  Nonlinear  Adaptive  Body  Axis  Washout  Circuit 
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FIGURE  C-22  - Time  History  For  Nonlinear  Adaptive  Body  Axis  Washout  Circuit  (Con't) 
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FIGURE  C-23  - Time  History  For  Nonlinear  Adaptive  Inertial  Axis  Washout  Circuit 
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FIGURE  C-23  - Time  History  For  Nonlinear  Adaptive  Inertial  Axis  Washout  Circuit  (Con't) 
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FIGURE  C-23  - Time  History  For  Nonlinear  Adaptive  Inertial  Axis  Washout  Circuit  (Con't) 
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FIGURE  C-24  - Time  History  For  Nonlinear  Adaptive  Washout  Circuit  With  Linear  Transformations 
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FIGURE  C-24  - Time  History  For  Nonlinear  Adaptive  Washout  Circuit  With  Linear  Transformations  (Con't) 
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Tho  axis  transformation  matrices  for  a link  three-degree-of-freedom  motion  base 
arc  given  as; 


where 


cese  0 -siNe 

SIN  ^ SING  cos<>  siN^ccee 

COS^SING  -SIN«  cos<>*  cose 


BODY 


BODY 


BODY 


INERTIAL 


INERTIAL 

^INERTIAL 


COS  e SIN p SIN  0 COS 0 sme 


COSO 


-SINO 


siNe  SINO  cose  cos  o cose 


where 


INERTIAL 


INERTIAL 


INERTIAL 


BODY 


BODY 


BODY 


The  transformation  from  pitch,  roll  and  heave  commands  (Oj^,  0^,  Z^)  to  right, 
left  and  rear  leg  commands  (L  RIGHT,  L LEFT,  L REAR)  Is  as  follows: 

Rear  Leg 


2.6025  -Z^  -AF^pSINej^ 


AX^G  = 7.5  10.3547  -SIN 


= X^Q  + (AF^p  + 2.3568)  (1  -COS  Gj^)  -1.875  SIN COSO 


L REAR 


AXp  + (2.6025  -(AF^p  + 2.3568)  SINGj^  -Z^ 


(1.875  COSe„  COS  0„)^  + (1.875  SIN  0 . -4.4775 

MM  M 
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Right  Leg 

= 

(4.81-  X^p)  (COSe^-1)  +AX^^-  3.075  SIN 

= 

2.6025  (-  (4.81  -■^F^p^  "^2 

TSP 

= 

= 

COS'^  (AZgp/TSP) 

TSP  SIN  (0j^^ -O 

^^R 

3.075  SIN  0 .,  +4.7344  (CCS  G -1) 

M M 

= 

TSP  COS  (Gj^^  -^)  -4.7344  SIN  0 ^ 3.075  (CCS  0 -1) 

^ RIGHT 

= 

>/  (AXp)^  + (AYp)2  + (AZp)2 

Left  Leg 

= 

3.075  SIN  0 . + 4.7344  (1  -COS  0 .,) 

M M 

^^L 

= 

TSP  COS  (Gj^  -4)  + 4.7344  SIN  0 ^^  + 3.075  (COS  0 -1) 

^ LEFT 

= 

\/  (AXp)^  + (AYj^)^  + (AZ^)^ 

A linearized  version  of  this  transformation,  accurate  to  better  than  5 percent  is 
as  follows; 

Rear  Leg 

^REAR  " ^2  +eM  ( F^p  + 2.3568) 

Right  Leg 

^lUGHT  " ^2  "®M  -t- 4.73440  ^ 

Left  Leg 

^LEFT  = 
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In  tlic  sc‘  transformations,  A F(^p  represents  the  distanee  in  feet  forward  to  the 
desired  eenter  of  pitch  from  the  pivot  point  attaching  the  simulator  to  the  floor. 

I’hese  equations  represent  the  leg  position  commands  in  feet  from  the  simulator  level 
position.  An  input  of  65.  844  volts  will  move  the  legs  one  linear  foot.  The  rear  leg  zero 
volt  position  is  also  the  zero  displacement  position.  A bias  of  30.864  volts  must  be 
applied  to  the  front  legs  to  bring  the  simulator  to  its  level  position. 
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APPENDIX  D 
ANALYSIS  TECHNIQUES 


The  statistical  tests  used  in  this  study  are  of  the  "sampled  permutation"  variety, 
an  extension  of  the  Fisher  Randomization  method.  The  sampled  permutation  technique 
was  chosen  because  it  requires  that  no  assumptions  be  made  concerning  the  underlying 
distribution  of  the  experimental  data. 

For  comparing  different  washout  types  and  the  use  of  washouts  in  different  axis 
systems,  (Phase  II  of  the  experiment)  the  statistic  used  was  an  F ratio.  Computation 
of  the  F ratios  for  row  differences  (washout  tjrpes)  and  column  differences  (use  of 
washouts  in  different  axis  systems)  was  carried  out  exactly  as  it  would  be  for  a standard 
two-way  analysis  of  variance  with  replication.  Ordinarily  these  F ratios  would  be 
compared  with  tabulated  values  for  the  F distribution  to  test  for  statistiqai  significance 
of  the  results.  Such  action,  however,  implies  an  assumption  that  the  e}q)crimental  data 
is  a sample  from  a population  with  a normal  distribution.  Rather  than  make  this  as- 
sumption, the  sampled  permutation  technique  generates  its  own  F distribution  from  the 
actual  experimental  data.  The  F ratios  for  rows  and  columns  are  then  compared  with 
the  newly  generated  F distribution  to  determine  the  statistical  significance  of  the  results. 

Table  D-I  contains  the  data  for  attitude  error  variances  from  the  second  phase  of 
the  ejqjeriment.  Analysis  of  this  data  will  serve  as  an  example  of  the  use  of  the  sam- 
pled permutation  technique. 

There  are  twelve  washout/axis  system  combinations.  Each  combination  has  four 
replications  - one  for  each  pilot  taking  part  in  this  phase  of  the  experiment.  Averages 
are  taken  over  the  four  pilots  for  each  of  the  washout/axis  system  combinations  to 
form  the  cell  means  of  table  D-n.  From  the  row  and  column  totals  in  table  D-II 
(washout  t3rpe  and  axis  system  of  application  respectively)  It  appears  that  the  attitude 
error  variances  are  roughly  equal  regardless  of  the  axis  system  used  while  the  vari- 
ances for  third  order  washouts  are  substantially  larger  than  for  the  other  washout 
types.  These  results,  however,  need  to  be  tested  for  statistical  significance.  The 
standard  two-way  ANOVA  table  for  this  data  is  shown  in  table  D-III.  At  this  point  the 
F distribution  must  be  generated  by  the  sampled  permutation  method.  The  data  from 
each  pilot  is  randomly  shuffled  and  a new  table  formed  from  the  shuffled  data.  F 
ratios  are  recalculated  and  checked  against  the  F ratios  from  the  original  data  arrange- 
ment. Once  again  the  data  from  each  pilot  is  shuffled  and  new  F ratios  calculated  and 
compared  with  the  original  F ratios.  This  procedure  is  repeated  many  tinies  and  a 
count  Is  made  of  how  many  times  the  actual  F ratios  exceed  those  generated  by  shuffling 
the  data.  The  size  of  the  data  sample  prohibits  examining  ail  permutations  of  the 
data  - thus  the  origin  of  the  term  "sampled  permutations".  The  implications  are  that 
the  F ratios  for  the  actual  experimental  results  will  exceed  the  F ratios  generated  by 
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TABLE  D-I 

ATTrrUDE  ERROR  VARIANCES  FOR  PHASE  II  OF  THE  EXPERIMENT 


Linear  Transformation 

Body  Axis 

Inertial  Axis 

First  Order 

.05869 

.03414 

.01370 

Washout 

.02152 

.03279 

. 03053 

.03058 

. 02223 

.05994 

. 04536 

.07591 

. 11633 

Third  Order 

. 12220 

.01938 

.01375 

Washout 

. 12968 

.07327 

. 13285 

.06578 

.04907 

.29300 

.06578 

. 07655 

. 04063 

Adaptive- 

.06357 

. 04340 

.01700 

Xonlirif.-ar 

.02664 

. 11704 

. 19751 

Washout 

.07545 

.01903 

. 0358S 

.05643 

.06711 

. 04235 

High  (jrder 

.02760 

.04384 

.01521 

Washout  with 

.01438 

.08191 

. 04337 

Tilt 

. 02982 

.03964 

.02213 

.05507 

.10487 

.04411 

TABLE  D-II 

CELL  MEANS  FOR  ATTITUDE  ERROR  VARIANCES  FROM 
PHASE  U OF  THE  EXPERIMENT 


Linear 

Transformation 

Body  Axis 

Inertial  Axis 

Row 

Totals 

First  Order  Washout 

. 15615 

. 16507 

. 22050 

.54172 

Third  Order  Washout 

.38713 

.21827 

.48023 

1.08563 

Nonlinear  Ad.aptive 

. 22209 

.24658 

. 29274 

.76141 

Washout 

High  Order  Washout 

. 12687 

. 27026 

.12482 

.52195 

with  Tilt 

Column  Totals 


89224 


90018 


1.11829 


2.91071 
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TABLE  D-ni 

ANOVA  TABLE  FOR  ATTITUDE  ERROR  VARIANCES 
FROM  SECOND  PHASE  OF  THE  EXPERIMENT 


I 


Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Squares 

F Ratio 

Rows 

.0171810 

3 

.0057270 

2.20 

Columns 

. 0020570 

2 

. 0010285 

.39 

Interaction 

.0114877 

6 

.0019146 

.74 

Subtotal 

.0307257 

11 

.0027932 

1.07 

Residual 

. 093764 

36 

. 00260446 

Total 

.12449 

47 

permuting  the  data  a large  percentage  of  the  time  if  the  results  are  due  to  actual  physi- 
cal differences  in  the  washout  algorithms  rather  than  just  ordinary  statistical  fluctua- 
tions. 


The  levels  of  significance 


® 


for  the  results  are  as  follows: 


Washout  Type  . 08 

Axis  System  of  Application  . 76 

The  implications  of  these  statistics  are  that  the  axis  system  of  application  for 
washouts  has  no  physical  effect  on  the  attitude  error  variances.  However,  there  is 
some  indication  that  the  poorer  performance  of  third  order  washout  filtei's  is  not  merely 
a statistical  fluctuation.  While  a significance  level  of  . 08  is  not  normally  considered 
convincing  (.  05  or  . 01  are  common  test  levels  for  significance)  the  author  feels  that 
this  is  a reportable  result  since  the  sampled  permutation  technique  is  not  considered  to 
be  as  sensitive  as  standard  analysis  techniques. 

For  Phase  in  of  this  experiment,  a two-way  analysis  of  variance  with  the  sampled 
permutation  technique  was  again  used.  Table  D-IV  shows  the  data  for  position  error 


Level  of  significance  here  refers  to  the  percentage  of  time  that  F ratios  generated 
by  randomly  permuting  the  data  exceed  the  F ratios  of  the  original  arrangement  of 
the  data. 
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TABLE  D-IV 

POSITION  ERROR  VARIANCE  FOR  PHASE  m OF  THE  EXPERIMENT 


First  Order 
Washout 

Adaptive  Nonlinear 
Washout 

High  Order  Washout 
with  Tilt 

Pilot  #1 

.4683 

1.6945 

.1655 

.2006 

.5818 

.2674 

.4857 

.3651 

.0562 

.3084 

.3432 

1.4548 

.2011 

1.5972 

. 1902 

Pilot  #2 

1.5214 

.6351 

.2412 

.4921 

.6001 

.8324 

.2214 

.7274 

4.3081 

.2929 

.2578 

.2193 

.5678 

.3053 

.3180 

Pilot  #3 

.3279 

.6875 

.4821 

.3080 

.1175 

.5671 

.0486 

.0684 

.2109 

.4269 

.1248 

.4025 

.2294 

.1214 

.0784 

variance  for  the  third  phase  of  the  e;q)erinient.  The  data  was  arranged  differently  for 
analysis  of  this  phase  of  the  e^qjeriment.  The  columns  now  refer  to  washout  types  and 
the  rows  contain  data  from  the  different  pilots.  Each  washout/pilot  combination  has 
five  replications.  The  ANOVA  table  is  set  up  as  before  and  the  sampled  permutation 
technique  is  applied  as  in  the  second  phase  of  the  e;q)eriment.  The  data  from  each  pilot 
is  shuffled  and  an  F distribution  is  generated  from  the  actual  data.  Gbce  again  the  im- 
plication is  that  the  observed  results  are  not  due  to  statistical  fluctuations  if  it  can  be 
shown  that  the  F ratios  calculated  from  the  original  arrangement  of  the  data  are  larger 
than  the  F ratios  generated  by  permuting  the  data. 

The  differences  in  performance  for  the  three  washout  filters  must  be  attributed  to 
statistical  fluctuations  rather  than  real  physical  differences  since  the  level  of  signifi- 
cance is  .75  in  this  case. 

individual  comparisons  of  washouts  for  the  third  phase  of  the  experiment  were  also 
made  with  a two-way  analysis  of  variance  with  sampled  permutation. 

A sampled  permutation  student  test  was  used  for  comparing  motion  versus  no 
motion. 
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